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Abstract 



io With local density approximation + Hubbard U (LDA+C7) calculations, we show that the fer- 
ii romagnetic (FM) insulating state observed in tensile-strained LaCo03 epitaxial thin films is most 



12 likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Com- 

13 pared with other FM states, including the intermediate-spin (IS) state (metallic within LDA+U), 

14 which consists of IS Co only, and the insulating IS/LS mixture state, the HS/LS state is the most 

15 favorable one. The FM order in HS/LS state is stabilized via the superexchange interactions be- 
i6 tween adjacent LS and HS Co. We also show that Co spin state can be identified by measuring the 
17 electric field gradient (EFG) at Co nucleus via nuclear magnetic resonance (NMR) spectroscopy. 



is PACS numbers: 75.50.Cc, 75.70.Ak, 76.60.Gv 



19 Perovskite-structure oxides have been proven a fertile area in condensed matter physics. 

20 They exhibit amazing properties, including ferroelectricity, ferromagnetism, colossal magne- 

21 toresistance (CMR), and multiferroics (simultaneous ferroelectricity and ferromagnetism), 

22 as a consequence of their spin, lattice, charge, and orbital degree of freedom. Advances in 

23 thin-film growth techniques have even brought more promising potentials for their future 

24 application, as their properties can be engineered via epitaxial strains. A few examples 

25 include strontium titanate (SrTiOs), ferroelectric in tensile-strained thin film while para- 

26 electric in bulk, 1 lanthanum titanate (LaTiOs), conducting in compressive-strained thin 

27 film while insulating in bulk, 2 and Europium titanate (EuTiOs), in which multiferroics in- 

28 duced by tensile strains has been observed. 3 As to lanthanum cobaltite (LaCoOs), a diamag- 

29 netic insulator in bulk at low temperatures (T < 35 K), a ferromagnetic (FM) insulating 

30 state has been observed in tensile-strained thin films, e.g. LaCoOs grown on SrTiOs or 

31 (LaA103)o.3(Sr 2 AlTa0 6 )o.7, at T < 85 K,-~— while the ferromagnetism induced by com- 

32 pressive strains, e.g. LaCoOs grown on LaAlOs, is not conclusive .— — « Two questions arise 

33 immediately: (1) Given that there are six 3d electrons in Co 3+ , which can thus have a to- 

34 tal electron spin S = 0, 1, or 2, refereed to as low-spin (LS), intermediate-spin (IS), and 

35 high-spin (HS) state, respectively, what is the spin state of Co in FM LaCoO"3 thin films, in 

36 contrast to the LS Co in diamagnetic bulk? (2) What leads to the formation of FM order in 

37 LaCoOs thin films? After all, FM insulators are rarely seen. So far, all first-principles cal- 

38 culations have only found FM metallic LaCoOs thin films with all Co ions in IS state,— ^^ 

39 a prediction clearly inconsistent with transport measurement s.- 

40 While finite Co spin induced by tensile strains in LaCoOs thin films has just started 
4i attracting attention, finite Co spin induced by thermal excitation in bulk LaCoOs has been 

42 a highly controversial issue for decades.— 1 ^ With LS Co at T < 35 K, bulk LaCoOs becomes 

43 a paramagnetic insulator with finite Co spin at about 90 K. Such a spin-state crossover in 

44 the temperature range of 35-90 K was first suggested to be a LS-HS crossover— ~— but was 

45 later suggested to be LS-IS based on a local density approximation + Hubbard U (LiDA+U) 

46 calculation.— Since then, both scenarios have received supports from various experimental 

47 and theoretical works, but a consensus is not yet achieved (see Ref. l23| for a brief review). 

48 A study regarding LaCoOs thin films may also help understanding LaCoOs bulk from a 

49 different perspective. In this paper, we investigate the Co spin state in tensile-strained thin 



50 films and the formation of FM order via a series of LDA+U calculations. While LDA+U has 

51 been frequently used to study cobaltites and Co spin state, the choice of Hubbard U can be 

52 an issue. It has been shown that under the same lattice parameter, the Hubbard U affects the 

53 total energy and the determination of ground state.— A well justified Hubbard U determined 

54 by first principles would thus be necessary for finding out the the actual ground state. In 

55 this paper, we compute the Hubbard U parameters of Co in all spin states self-consistently 

56 with a linear response approach.— ~— This method has successfully found the ground state 

57 of iron-bearing magnesium silicate (MgSiOs) perovskite at a wide range of pressure.— Both 

58 the plane-wave pseudopotential (PWPP) method 2 - implemented in Quantum ESPRESSO 

59 codes 28 and the augmented plane wave + local orbitals (APW+lo) method 29 implemented 
eo in WIEN2k codes 3 -^ are used. As shall be pointed out later, the orbital occupancies of Co in 
ei thin films are different from those in bulk, due to their different symmetries. We therefore 

62 compute the electric field gradient (EFG) tensor at Co nucleus, V zz , with WIEN2k, to see 

63 whether the Co spin state in thin films can be identified via EFG, as demonstrated in bulk. 23 

64 The pseudocubic lattice parameter of bulk LaCo0 3 (R3c symmetry) is about 3.81 A 

65 at T ~ 5 K.— »2S To model tensile-strained LaCoOs thin films via bulk calculations, we 
ee constrain the in-plane pseudocubic lattice parameters a pc and b pc of the hypothetical bulk 
67 to 3.899 A (the lattice constant of cubic SrTiOs at low temperatures),- 3 set a = (3 = 7 = 
6B 90°, and optimize the out-of-plane pseudocubic lattice parameter (c pc ). Due to the lack of 

69 accurate information regarding C0O6 octahedral rotation in thin-film LaCoOs in the low- 

70 temperature FM phase, we consider two extreme cases shown in Fig. [TJ (a) cube-on-cube, 

71 namely, no CoOg octahedral rotation degree of freedom, and (b) full CoOg rotation degree 

72 of freedom subject to the above-mentioned constraints. We also consider several magnetic 

73 configurations shown in Fig. [2j (a) all Co ions in LS state, (b) all Co ions have the same 

74 magnetic moment aligned in FM order, and (c) a mixture state with LS Co surrounded by 

75 magnetic Co (and vice versa) aligned in FM order. The configuration shown in Fig. 2(c) is 

76 a legitimate postulate, as the observed magnetization in LaCoOs thin films rarely exceeds 

77 0.85 hb/Co.—~— For the configuration in Fig. [2(b), a convergent wave function for HS state 
7B cannot be obtained; only IS state can be found. For the configuration in Fig. [2(c), both 
79 HS/LS and IS/LS mixture state can be obtained. The self-consistent Hubbard U parameter 
so (U sc ) of LS and IS Co are 7.0 eV, while the HS/LS state has f/ s ( c HS) = 5.4 and u£ S) = 7.2 



si eV.— The dependence of U sc on c pc is negligible. To demonstrate how the choice of Hubbard 

82 U can affect the determination of ground state, we also present the result obtained using 

83 a constant U = 7 eV for all Co in our PWPP calculations. In tensile-strained LaCo03 

84 thin films, CoOg octahedra possess tetragonal symmetry, namely, longer Co-0 distance on 

85 the xy plane, regardless of CoOg rotation, as will be discussed in the next paragraph. In 

86 tetragonal symmetry (Ai/i), the spin-down electron of HS Co occupies d xy orbital, and the 

87 spin-down electrons of IS Co occupy d xz and d yz orbitals, as shown in Fig. [2](d) . Such orbital 
as occupancies are very different from those in bulk LaCoOa, in which Co06 octahedra have 

89 trigonal symmetry (D^d), with the [111] direction being the high-symmetry axis. In bulk 

90 LaCoOa, the spin-down electron of HS Co occupies the 0^2-like orbital oriented along the 

91 [111] direction, and the IS Co spin-down electrons occupy the doublet with 3-fold rotation 

92 symmetry about the [111] direction.— 

93 The optimized out-of-plane pseudocubic lattice parameter (c pc ) of each FM state and 

94 associated relative energy (AE) and band gap (E gap ) are listed in Tables I-II. Regardless 

95 of CoOq rotation, the HS/LS mixture state (with U sc ) is the most stable FM state given 

96 by the PWPP method (Table l)M While the choice of U = 7.0 eV makes HS/LS state 

97 less favorable in PWPP calculations, APW+lo calculations still find HS/LS the most stable 

98 FM state (Table II). Both PWPP and APW+lo methods open an energy gap for HS/LS 

99 state, consistent with transport measurements. 6 Also, the presence of HS Co is consistent 
ioo with recent x-ray magnetic circular dichroism (XMCD) and x-ray absorption spectroscopy 
ioi (XAS) spectra.- 1 ^ In contrast, the IS state is never the most favorable FM state, regardless 

102 of the computation method and CoOg rotation. Its partially filled bands formed by partially 

103 occupied d z 2 and d x 2_ y 2 orbitals in IS Co lead to a nonzero density of state at the Fermi 

104 level. When the IS Co concentration is reduced to 50%, an energy gap is opened in APW+lo 

105 calculations (Table II), while PWPP calculations still give a conducting IS/LS state (Table 
we I). Such a difference is likely to result from the way that the Hubbard U is applied. In 
107 PWPP, the Hubbard U is applied to the projection of the total wave function onto Co 3d 
io8 orbitals; 24 in APW+lo, the Hubbard U is directly applied to the 3d orbitals within the 
109 muffin-tin radius of Co (1.9 bohr). Insulating or not, the IS/LS state is highly unlikely; 
no its energy is even higher than that of IS state. One more thing worthy of mention is that 
in the Co-0 distance (rfco-o) an d Co-O-Co angle obtained in our calculation are different 



112 from those estimated in Ref. |5|, where a constant dc -o = 1-93 A is assumed (regardless 

n3 of compressive or tensile strains), and an Co-O-Co angle of 176° is estimated for LaCoOs 

H4 grown on SrTi03. In our PWPP calculation with U sc and CoOg rotation, the HS/LS state 

us has dco(HS)-o — 2.015 and 1.872 A, g?co(ls)-o = 1.922 and 1.886 A, and Co-O-Co angles of 

n6 163.8° and 156.6°, on the xy plane and along the z axis, respectively. Even for the IS state 

H7 suggested by Ref. |5|, we have (Zco(iS)-O — 1-973 and 1.939 A, and Co-O-Co angles of 162.3° 

us and 154.9°, on the xy plane and along the z axis, respectively. 

n9 Other than the total energy, structural properties can be a useful criterion to determine 

120 which FM state favors tensile strains (c pc /a pc < 1). Starting with the structures listed in 

121 Table I with CoOg rotation, we perform full structural optimization (at constrained volume) 

122 via variable cell-shape damped molecular dynamics. 3 - All lattice parameters, including a, 

123 /3, and 7, are optimized, so the final structures only experience hydrostatic pressures. With 

124 a, (3, and 7 slightly deviated from 90°, the IS state has c pc /a pc > 1, while all other states 

125 remain c pc /a pc < 1 (but IS/LS state still has a c pc /a pc larger than that of HS/LS state), as 

126 shown in Table III. The larger c pc /a pc ratio associated with IS Co is a direct consequence of 

127 its occupied d xz and d yz orbitals (by spin-down electrons), which elongate the Co-0 distance 
i2B along the z-direction. In contrast, the fully optimized HS/LS state has c pc /a pc = 0.969, in 

129 great agreement with c pc /a pc = 0.967 observed in experiments- 

130 A significant part of cobalt-spin controversy arises from the difficulty in directly measuring 

131 the total electron spin of Co. Such difficulty, also appearing in other spin systems, can be 

132 resolved by comparing the calculated and measured EFGs.— -2SM. So far, insulating FM state 

133 has been observed in LaCoOa thin films with a pc ranging from 3.84 to 3.90 A.-<£ In these 

134 thin films, the magnetic Co concentration and Co-0 distance may be different, which can 

135 lead to slightly different EFG for Co in the same spin state. To find out possible upper 

136 and lower limits of HS and IS Co EFG, we compute them in two extreme cases: (1) thin 

137 films with a pc = 3.899 A and 50% of magnetic Co, namely, the HS/LS and IS/LS states 
i3B listed in Tables I and II, and (2) single isolated HS or IS Co in an array of LS Co in a fully 

139 relaxed structure with a pc ~ 3.81 A, where the orbital occupancies of isolated HS and IS Co 

140 are maintained in tetragonal symmetry [Fig. |2{d)]. In these APW+lo calculations, IS Co 

141 does not lead to a metallic state, in contrast to bulk LaCoOs {D^d symmetry).— Different 

142 choices of Hubbard U have been adopted as well (5 eV, 7 eV, and U sc ). The results of all 



143 these calculations show that the EFG mainly depends on the spin state: 14.7 < V zz /(10 21 

144 V/m 2 ) < 19.9 and -14.6 < \4 IS) /(10 21 V/m 2 ) < -8.0. The quadrupole frequency, v Q = 

145 3eQ\V xx \/2I (21 - 1) h, can thus be easily predicted, with Q = 0.42 x 10~ 28 m 2 and I = 7/2 
we for 59 Co nucleus. Based on the range of V zz and V zz , we conclude that in insulating 
147 LaCoOs thin films, a measured vq via nuclear magnetic resonance (NMR) spectroscopy 
us within ~ 8.2 ± 2.4 MHz can be a strong evidence for IS Co, and a measured vq within 

149 ~ 12.6 ± 1.9 MHz should indicate HS Co. 

150 Analysis of electronic structures can help developing a physical understanding for the FM 

151 order in HS/LS state, whose projected density of states (PDOS) are shown in Fig. |3j The 

152 case with U = 7 eV and no CoOg rotation is presented, as the main features in PDOS are 

153 not sensitive to the choice of U and CoOg rotation. Extracted from Fig. [3]^a), both HS and 

154 LS Co have nonzero magnetic moment: 2.97 and 0.56 p#, respectively. The FM order is 

155 established via the superexchange interaction between HS and LS Co, as described by the 

156 Goodenough-Kanamori rule,— >^~— which states that the superexchange interaction between 

157 two cations (with or without a shared anion) is ferromagnetic if the electron transfer is from 

158 a filled to a half-filled orbital or from a half-filled to an empty orbital. Indeed, for the HS/LS 

159 state, electrons transfer from the filled d xz and d yz orbitals of LS Co to the half-filled d xz 
wo and d yz orbitals of HS Co via the oxygen in between, and also from the half-filled e g (d z 2 
lei and d x 2_ y 2) orbitals of HS Co to the empty e g orbitals of LS Co, as depicted in the inset of 

162 Fig. [3]^b). The PDOS shown in Fig. EJ^b) confirms this model: the finite spin- up e g electrons 

163 localized at the LS Co site (transferred from the HS Co site) and the finite spin-down d xz 

164 and d yz electrons localized at the HS Co site (transferred from the LS Co site). Such electron 

165 transfers have been also described via a configuration fluctuation model,— ^ which further 
we suggests that the interchange of spin states (without net transfer of charge) led by electron 
167 transfers stabilizes the FM order in this HS/LS states. 

lea The above-mentioned superexchange interaction can be visualized via electron spin den- 

169 sity s(r) = p-|-(r) — p±(r), where Pt( r ) an d Pi( r ) are spin-up and spin-down electron density, 

170 respectively. Figure Hl^a) shows s(r) corresponding to the configuration with all HS Co 

171 magnetic moments aligned (same as the configuration in Fig. [3]). The nonzero magnetic 

172 moments localized at LS Co sites (with e g character), aligned with the HS Co magnetic 

173 moments, are consistent with the PDOS shown in Fig. [3j When the magnetic moment of 



174 one HS Co in a 40-atom supercell is nipped [Fig. 11(b)], the alignment of magnetic moments 

175 is altered, and so is the condition that allows configuration fluctuation [inset of Fig. 13(b)]. 

176 The spin density at surrounding LS Co sites is thus significantly affected. One flipped HS 

177 Co spin (in a 40-atom cell) increases the total energy by 195 meV/supercell. Flipping one 
i7B more HS Co spin, so the total magnetization per supercell becomes zero, further increases 
179 the total energy by 78 meV/supercell. With CoOe rotation, the energy increases associated 
i8o with one and two flipped HS Co spins are 96 and 34 meV/supercell, respectively. These 
i8i results indicates that the magnetic moment of HS Co in the HS/LS state should align at 

182 low temperatures. 

183 While our calculations have shown that HS/LS state is the most favorable state among 

184 the ferromagnetic states being considered, magnetic state in actual LaCo03 thin films 

185 can be more complicated. The magnetization observed in experiments rarely exceeds 0.85 

186 fiB /Co,— — smaller than that of the HS/LS mixture (2 //b/Co). Such a magnetization sug- 

187 gests that the HS Co population should be smaller than 50%. In fact, XAS spectra combined 

188 with atomic multiplet calculations have suggested that LaCoOs thin film on SrTiOs consists 

189 of about 64% of LS Co and 36% of HS Co.— Given that the FM order is achieved via the 

190 superexchange interaction within the HS-LS-HS Co configuration shown in Fig. E(c), one 

191 can thus expect ferromagnetic HS/LS domain and nonmagnetic LS domain coexist in tensile 

192 strained thin films, as observed using magnetic force microscopy (MFM). 7 Also, since HS/LS 

193 state favors larger in-plane lattice parameter, thin films with larger in-plane lattice parame- 

194 ters can be expected to have larger HS/LS domain, and thus larger magnetization, consistent 

195 with the increase of magnetization with lattice parameter observed in experiment .— 

we In summary, we use LDA+U calculations to show that the ferromagnetic insulating state 

197 in tensile-strained LaCoOa thin films is most likely a mixture of HS and LS Co. Among 

198 all the ferromagnetic states studied in this paper (HS/LS, IS/LS, and IS), the insulating 

199 HS/LS mixture state is the most favorable one, energetically and structurally. Its FM order 

200 is established via the superexchange interaction between LS and HS Co. We also show that 

201 cobalt spin states in LaCoOa thin films could be identified via NMR spectroscopy. 
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TABLE I. Optimized out-of-plane pseudocubic lattice parameter (c pc ) and associated relative en- 
ergy (AE) and energy gap (E gSu p) of each FM state in tensile-strained LaCo03 thin film (PWPP 
method). 

No CoOg rotation 





Cpc (A) 


AE 


(eV/f.u.) 


E ga ,p (eV) 


LS 


3.865 




0.35 


0.54 


IS 


3.785 




0.20 


metal 


IS/LS 


3.720 




0.35 


metal 


HS/LS (U sc ) 


3.685 




0.00 


0.92 





Full CoOg rotation 


Cpc (A) 


AE 


(eV/f.u.) 


£ga P (eV) 


3.660 




0.32 


1.24 


3.785 




0.19 


metal 


3.680 




0.29 


metal 


3.680 




0.00 


0.90 



HS/LS (U = 7 eV) 3.700 



0.29 



1.12 



3.695 



0.29 



0.90 



TABLE II. Optimized c pc and associated AE and E gap of each FM state (APW+lo method, with 
CoOg rotation). 







Cpc (A) 


AE 


' (eV/f.u.) 


£gap (eV) 


LS 




3.660 




0.37 


1.72 


IS 




3.741 




0.18 


metal 


IS/LS 




3.672 




0.29 


0.59 


HS/LS (U = 7 


eV) 


3.686 







1.52 



TABLE III. Fully optimized pseudocubic lattice parameters of each FM state at the volume as in 
Table I (with CoOg rotation). 







%>C7 "pC (A) 


Cpc (A) 


Cpc/Opc 


IS 




3.848 


3.893 


1.012 


IS/LS 




3.847 


3.778 


0.982 


HS/LS (Use 


) 


3.863 


3.745 


0.969 


HS/LS (U = 7 


eV) 


3.865 


3.757 


0.972 



10 




FIG. 1. (color online). Possible atomic structures of LaCoO"3 thin film (La is not shown) subject 
to constrained in-plane lattice parameters, (a) Cube on cube, no C0O6 octahedral rotation; (b) 
full octahedral rotation degree of freedom. 





FIG. 2. (color online), (a)-(c) Possible magnetic configurations in LaCo03 epitaxial thin film (La 
is not shown). The arrows denote for nonzero magnetic moments, either IS or HS. (a) LS state; (b) 
HS or IS state in FM order; (c) HS/LS or IS/LS mixture state in FM order, (d) The 3d orbitals 
occupied by the spin-down electrons of HS and IS Co. 
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Energy (eV) 



FIG. 3. (color online). Projected density of states of ferromagnetic HS/LS state (no CoOg rotation, 
£7 = 7 eV) onto (a) each atomic site, and (b) some of the Co 3d orbitals. The inset in (b) shows 
the electron transfer between HS and LS Co. 




FIG. 4. (color online). Spin density, s(r), of HS/LS state (no Co06 rotation, U = 7 eV) with 
(a) all HS Co magnetic moments aligned, and (b) one HS Co magnetic moment flipped downward 
(indicated by arrow). The isosurface values are 0.02 (red) and —0.02 (blue). 
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